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ABSTRACT 


y)  This  paper  Is  the  first  one  in  the  series  of  three  which  are 
addressing  in  detail  the  properties  of  the  three  basic  versions  of  the 
finite  element  method  in  the  one  dimensional  setting.  The  main  emphasis 
is  placed  on  the  analysis  when  the  (exact)  solution  has  singularity  of 
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1 .  INTRODUCTION 

In  general,  the  h-verslon  of  the  Finite  Element  Method  is  the 
standard  version  when  the  degree  of  elements  is  fixed  and  the  convergence 
is  achieved  by  the  refinement  of  "he  mesh.  The  p-verslon  fixes  the  mesh 
and  the  convergence  is  obtained  by  the  increase  of  the  degree  of 
elements.  The  h-p  version  simultaneously  refines  the  mesh  and  increases 
the  degree  of  elements.  In  recent  years,  the  p  and  h-p  versions  of  the 
finite  element  method  attracted  large  interest  both  in  theory  and 
computational  practice.  The  first  commercial  program  PROBE  based  on  the 
p-version  become  available.  This  paper  is  the  first  in  the  series  of 
three  which  will  analyze  in  detail  the  properties  of  the  h,  p  and  h-p 
versions  for  solving  the  one  dimensional  problem,  the  solution  of  which 
has  singularity  of  xa-type.  The  first  part  analyzes  the  p-version,  the 
second  the  h-version  and  general  h-p  versions  and  the  final  third  one 
addresses  the  problem  of  the  adaptive  h-p  version. 

In  this  paper  we  will  not  discuss  the  two  dimensional  case;  neverthe¬ 
less  the  detailed  results  in  one  dimension  are  serving  as  guidelines  to 
the  two  dimensional  theory.  For  the  analyses  of  the  p  and  h-p  finite 
element  version  for  the  two  dimensional  problem  we  refer  to  [1].  For  the 
computational  Implementation  and  engineering  aspects  of  the  p  and  h-p 
versions  we  refer  to  [2]. 

The  p-version  was  studied  theoretically  first  in  [3].  The  h-p  version 
was  addressed  in  [4],  and  for  detailed  analysis  of  the  p-verslon  in  three 
dimension  we  refer  to  [5]  [6].  For  additional  theoretical  aspect  of  the 
p-version  we  refer  also  to  [7].  The  p-verslon  is  in  some  sense  related  to 
the  theory  of  spectral  method.  For  these  results  we  refer  especially  to 
[8]  [91  [10J. 
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In  this  paper  we  will  consider  the  most  simple  model  problem 
(1.1)  -u"  -  f 

u(0)  -  u(l)  -  0 

with  the  solution 


(1.2)  uQ(x)  -  (x-?)“  -  (1-0°  x  -  (-£)“(l-x) 


where 


(*-£)“ 


(x-O 


If  x  >  5 


if  x  <  e. 


We  will  be  Interested  In  the  accuracy  of  the  finite  element  method 
measured  In  the  energy  norm  lul£  -  /q  (u')  dx.  We  have  to  assume  that 
a  > 1^  to  get  finite  energy  of  the  solution.  The  energy  norm  Is  In  our 
case  obviously  equivalent  to  the  H^-norm  of  the  standard  Sobolev  space 
H1. 

The  xa-type  singularity  of  the  solution  Is  an  analogue  of  the  singu¬ 
larity  of  the  solution  of  the  two  dimensional  boundary  value  problems  for 
elliptic  partial  differential  equation  occurlng  when  the  domain  has 
corners.  Solution  of  this  type  belongs  to  low  order  Sobolev  spaces. 

This,  in  general,  leads  to  low  rate  of  convergence  of  the  finite  element 
method  with  quaslunlform  meshes.  On  the  other  hand,  taking  Into  account 
the  special  structure  of  u(x),  one  can  achieve  high  (exponential)  rate 
of  convergence  by  proper  design  of  the  mesh  and  degrees  of  elements. 
Obviously,  the  results  related  to  (1.1)  can  be  generalised  to  the  general 
case  of  two  point  boundary  value  problem. 

For  our  model,  the  finite  element  method  with  C°  elements  gives  exact 
solution  in  the  nodal  points  and  the  analysis  of  the  error  of  the  finite 
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element  method  reduces  to  the  analysis  of  the  best  Lj-approximation  of 
u'  by  piecewise  polynomials. 

Most  of  our  results  In  this  paper  will  provide  both,  the  upper  and 
lower  bounds  of  the  error  of  the  finite  element  solution.  In  what  it 
follows  the  following  notation  will  be  used.  By  A  s  B  we  mean  that  A 
asymptotically  equals  B  as  some  parameter  tend  to  a  certain  limit.  By 
A  •  I  we  mean  that  A  Is  equivalent  to  B,  that  is,  there  exists  an 
equivalency  constant  C  >  0  (depending  on  some  parameters  to  be 
Indicated)  such  that 

CB  <  A  <  C_1B. 

The  paper  Is  organized  as  follows:  In  Sections  2-4  we  study  the  prop¬ 
erties  of  the  Legendre  expansion  of  the  function  (x-£)+-  Section  5 
addresses  the  Lengendre  expansion  of  more  general  functions.  Section  6 
discusses  the  performance  of  the  p-version,  and  Section  7  deals  with 
numerical  computations. 


2.  THE  LEGENDRE  EXPANSION  OF  (x~0+. 
Let 


„  ,  1  d  r/  2  ,  Nn, 

p  (*)  “  — - -  l(x  -1)  ] 

q  « n  *  ,  n 

2  n!  dx 

be  the  Legendre  polynomials  which  form  an  orthogonal  basis  of  L2[“l»l]< 
For  the  properties  of  Legendre  polynomials  see  [11],  [12],  [13],  [14]. 
Ue  denote 

(2.1)  u?(x)  -  (x-C)“. 

Suppose  that  its  Legendre  expansion  is 

00 

u  (x)  ~  l  a  P  (x). 


n**0 


n  n 


(2.2) 


The  following  theorem  gives  the  expression  for  the  coefficients  of  the 


t 


(2.8)  •  (x)  -  F( n-a, -a-  ^  n+  y;x) 

n,a  l  i 


oo  (  -a-  — )  (  — a+  — )  •  •  •  (  -a-  +k— 1  } 

1  +  y  (n~a)(n~ct+l  )•  •  *  (n-crHc-1 )  v  2M  2;  V  2  Ki;  k 

k-1  (n+  |)(n+  |)-*-(n+  |  +  k-1) 


00  k 


k! 


1  +  1  t  ~  ©££-)]  »k. 

k-1  j-1  ^  '2  * 


Proof.  We  have 


(2.9) 


2n+l 


a  -  -=■  ~  /  u  (x)P  (x)dx. 

n  2  1,  ^  n 


-1 


In  the  case  C  -  -1 


2n+l 


/  (l+t)aP(t)dt 


/  (l+t)a  [(t2-l)n](n)  dt. 

2  -1  2nn! 


I  I” 


Since  a  >  -1,  for  k  *  l,2,...,n  we  have 


llm  (l+t)oHc+l[(t2-l)n](n_k)  -  0, 

t+±l 


Integration  by  parts  yields 


a  -  (-1)"  a(a-l )» • » (a-iH-1 )  /  (l+t)a(t-l)n  dt 

n  2n+In!  -1 


2n+l  (-1  )2n_tt(a~l )  •  •  *  (a-n+1 )  }  n+t>«+n  d 
2rtfl  (a+l)(a+2)...(a+n)  ^  az 


it* 


- M.. 
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m  2n+l  a(ct-l )« « « (g-n+1 )  .ar+n+1 

2n+l  (a+l)(a+2)»*«(a+n+l) 


and  (2.3)  follows. 

For  the  case  -1  <  E  <  1  we  have 

(2-10)  an  -  (t-€)a  Pn(t)dt. 


For  Re  (x  >  0  and  0  <  y  <  1  one  has 


(2.11)  /  (I-.)**-1  P„(l-T»)&  -  p^’ci-,) 


(see  [14],  p  833).  Setting  t  *  1  -  (l-£)y  we  obtain  from  (2.10): 


an  "  ^2±L/i  [Cl-C)(l-y)]a  Pn(l-(l-5)y)[-(l-F)]dy 


-  «-o“+1 /  ci-,)^1)-1  P  (i-u-t),)d, 

0  n 


_  ?n+l  ,,_,,a+l  T(a+l)n!  „(a— 1 ,-a-l ) 
2  U  F)  r(a+n+2 )  Pn 


-  2trn  (l-F^ni  (a+l,-a-l) 

2  (a+l)(a+2)»»«(a+rt+-l)  n 


This  proves  (2.4). 

Finally,  for  F  <  -1  integration  by  parts  gives 


an  "  /  (‘“5)°  Pn(t)dt 


/  (t+UI)a— i— 
2  -1  2n.nl 


[(t2-l)n](n)dt 


9 


(-1)2  *“±L  fl  (t+k  | )«-n(t2_1  }n 


2n.n! 


dt 


Setting  t  *  cos  6,  we  obtain 


(-l)n  /  (t+lE; | )ac-n(t2-i)n  dt 
-1 


it  .  2n+l  „ 

/  - L_de. 

0  (l^l+cos  e)n  a 


by  [14],  p.  384  we  have 


it  .  2u-l 

j  _ sin  _  x 

0  (l+2a  cos  xfa^) 


2  v  dx  "  V2  )F(v  ,v-p+ V2  ;p+ V2  ;a2) 


with  Re  (i  >  0,  lal  <  1,  and  F(a,p;v;x)  being  the  Gauss  hyper geometric 
function.  Let 


kl  +  a2-i 


thus 


(0  <  r  <  1) 


|F,| 


1+r 
2r  ' 


Setting  H  *  n  +  1,  (p  >  0)  and  v  *  n  -  a,  we  have 


u  .  2n+l„ 
J  sin  8 


0  (l^l+cos  9) 


„  n  ,  2(n+l  )-l  „ 

.  (2„—/  — sis - J-Lde 

0  (1+2 r  cos  9+r  ) 


_  _sn-a  r(tt+l)/it  13? 

(2r)  , - 3T  F(n~«.-a  -  7; "  +  ~;r) 


rfn  + 


(2r) 


n-a  2 11+1  n! 


IT*  „<*> 


(2n+l)i!  n,a 


Wlth  $n,a  defined  In  (2.8).  (2.6)  now  follows  easily. 
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Remark  1.  For  5  “  -1,  the  expansion 


(2.12)  (l+t)°  -  2af— +  y  (  ccit+12(  2n+lJ  p  ,  >i 

^  ;  z  ta+1  +  L  ( a+1 ) ( a+2 ) •  • » ( a+n+1 ) 


Is  well  known  as  Neumann-Stiel jes  series  (see  [12],  pp.  240-244). 


Remark  2.  Since 


(2.13) 


ot(a+l )  •  •  *  (cr+n)  -  r(^-^ 

”(a) 


the  above  results  can  also  be  written  in  the  following  form: 
1)  if  £  -  -1,  a  >  -1,  then 


(2.14) 


/■  _i  \ n  T ( a+1 ) T ( n-a )  _na 

(  U  n-aTrTiSa+2)  (2n+1)  2 


.n-1  fr(l+a)1  sin  na  r(n-a) 


f ( a+n+2 ) 


2)  if  -1  <  £  <  1,  a  >  -1,  then 


(2.15) 


r(i+q)(i-n  r( 

r (a+n+2) 


(n  +  i>pi°t+1’~a~1)<?) 


3)  if  C  >  -1,  then 


(2.16) 


r(a+i)r(y  ) 

a  -  (-1)°  - — - — r-  (2r)n  a  ♦  (r2) 

T(-a)2  T(n  +  j)  n,a 


(_l)n“l  r(l+a)sln  Tta  r(n-a)  fn-a  .  2. 

2°V7  r(n  +  j)  C  n,a 


Remark  3.  In  the  special  case  £  ■  0,  we  have  (a  >  -1) 


(2.17)  a 


2n 


.  1.  (-l)ntl2-(a+l)2][32-(a-H)21»»«[(2n-l)2-(g-H)2] 
2>  (a+1 )(a+2 ) • • • (a+2n+l) 


2  1  rn+rr)  22nr(n+l  +  f)r(n  -  f) 

(-l)2(2n  +  b  - Z&& - FT?-™ - ~ 

2  (l+f)r(-f>  r(2n+ct+2 ) 


(2.18)  a2n+1 


(2n  +  I,  (-Dni22-(a+l)2l[42-(g4l)2J...[(2n)2-(a-K)2I 
'  2'  (a+2)(a+3)»»*(a+2n+2) 


(-l)2(2n  +  h 


-xa+°J 


2  r(-—)rd  -  ^±1) 


22n+1r(n+l+  ~)r<n  -  ~) 
r  ( a+2  n+3  )  ' 


Proof .  In  fact,  for  the  Jacobi  polynomial  P^a+*’  a  ^(t)  the  following 


recursion  formula  (see,  for  example  [11],  p.  71)  holds: 


(2.19)  P(a+1*~a  !)(t)  ,  i 

p(a+l,-a-l)(t)  „  i  +  a  +  t 


(a+l,-a-l),  .  _  2n+l  -,(g+l,-g-l),  .  .  n2-(q+l)2  p(tt+l,-a-l),  .. 

n+1  n+1  n  n(n+l)  n-1  (C) 


In  particular, 


(2.20)  p(a+l,-a-l)(0)  -  1 


P 


(a+1 ,-a~l ) 
1 


(0) 


1  +  a 


(a+l,-a-l) 

n+1 


(0) 


n2-(a+l)2  p(a+l,-a-l) 
n(n+l)  n-1 


Equations  (2.17)  and  (2.18)  are  obtained  by  discussing  the  cases  for  n 
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even  and  n  odd  reap. 


Remark  4.  Equations  (2.17)  and  (2.18)  can  be  combined  into  one  formula: 


(2.21) 


[y]  r(Htt)4  (a)  ,  2nr(-^+ Dr(-~) 

(_1)  IT  (n  2)  r  ( n+ct+2 ) 


with 


cos 


2  ’ 


<t>  (a) 


.  ita 

I  sin  , 


If  n  Is  odd 


If  n  Is  even 


and  [x]  being  the  largest  Integer  which  is  not  greater  than  x. 


3.  ESTIMATES  AND  ASYMPTOTIC  BEHAVIOUR  OF  THE  COEFFICIENTS  OF  THE  LEGENDRE 
EXPANSION 

In  this  section  we  will  obtain  the  asymptotic  formulae  and  some 
estimates  for  the  Legendre  expansion.  First,  one  can  easily  prove  the 
following  lemma  by  using  Stirling's  formula. 


Lemma  1 . 


(3.1) 


r(_n+a)  _  _1 —  „  0(i)) 
r ( n+B )  B-a  V  V' 


(n  ♦  «>) 


where  0{— )  depends  on  a,  6. 

We  now  prove  the  following  theorem  which  is  concerned  with  the 
asymptotic  behaviour  of  the  coefficients  obtained  in  Section  2. 


Theorem  2.  Let  an  be  the  coefficients  of  the  Legendre  expansion  of 
(x-£)+  and  a  > 


Then 


13 


1) 

(3.2) 

with 

2) 

(3.3) 

with 

3) 

(3.4) 

where 

n 

Proof . 
for  £ 

For  - 


if  l  -  -1,  then 


a 

n 


(-1) 


n-1  C0(a) 
2a+l 


(i  +  0(^)1 


C0(a) 


2g+1r(l-ha)^sln  ic a 


if  -1  <  F,  <  1,  then 


a+  V; 


an  "  /|r(a+l)f^~)  %osf(n  +  ^)6-(a  +  |)£]  +  0(1)1 


0  =*  arc  cos  Fa,  and  0(1)  holding  uniformly  for  Itl  <  1  -  e,  e  >  0, 
if  ^  <  -1,  then 


an  “  (-1) 


n-1  Cl(a) 


2  “  +  2  .  _,1 


a  + 


J  [  (1-r  )  +  0(^)1  •  r 


n-a 


c  (a)  ■  r,(ita)siiL^>  r  , - I -  and  o  m  — — — 1  >  0, 

2V,  ki+/^r  «+! 

holds  uniformly  with  respect  to  r  €  (0,1). 


By  Lemma  1  and  (2.14),  (3.2)  follows  immediately: 

-  -1 

*„  '  '■1>rV'fWW('*!) 

n 

1<£<1,  we  use  the  following  asymptotic  formula  for  Jacobi 
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1 


I  I'' 


polynomial  (see  [11]  p.  196). 


(3.5)  P(ct,P)(cos  9) 
x\ 


cos(N9+y) 


_  “  +  7  p  +  — 

/nn(sin  -j)  (cos  j) 


Y  +  Ofn  ) 


where 


n  +  (a+0+l)/2 


•<“  +  T>  f 


0  <  0  <  it. 

This  formula  holds  for  arbitrary  real  a  and  0  and  holds  uniformly  for 
0  (  [e,it-e],  £  >  0.  Thus,  in  particular,  setting  0  *  arc  cos  K,  we  obtain 

/_+!  _„_n  cosf  (n  +  i)Q-(a  +  i-)f] 

p(a+l,-o-l)(  .  - 2  - 2_2 -  +  ofn  2] 

n  ,  j  i 

_  Qa+T  a  ~a  ~  T 

/itn(sin  y)  (cos  j) 


sin  0)  ^  cos[  (n  +  y)9-(a  +  y)5-]  -=• 


m 


/nn  (1-cos  0) 


<r+l 


2  2J+0fn2). 


Since  (l-£)a+^  *  (l~cos  0)a+*,  by  (2.15)  we  get 


_  r(l+a)(l-Oa+1r(n+l)  ,  .  u  _(a+l  ,-a-l )  ,_* 
r ( n+a+2 ) -  (n  +  V  Pn  (° 


.  r(l+a)(l~coa  ^(np^))--^  (.sln_9) 
n  /it 


a  +  i 


1  i  ,  - 

2  cos[ (n  +  ±-)0-(a  +  f)f] 


(1-cos  9) 


0+1 


I 


r~ 
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a+l 

(l  r(1+a)(slS_i)  2  {cos[(n  +  I)  -  (a  +  h  +  o(i)}. 

✓  7t  n  2  2  2  n 


This  holds  uniformly  In  9  €  [e,it-e]. 

For  F  <  -1,  (2.16)  and  Lemma  1  gives  (uniformly  In  F) 


a  -  (-l)""1  C.(g)  — a  (r2) (1+0(1)). 

n  11  n,a  n 

a  +  y 

n 


It  remains  to  analyse  $  (r  ).  The  series 
'  n.a 


CO  CD 


*  +  T,  2k 


♦„  -  1  +  MTT  2)  r 

,a  k-1  3-1  t*j  +  i  k 


converges  uniformly  In  r  f  (0,1)  since  Its  general  term  is  dominated  by 


a  +  y 

(  2i 


i  r>  | 

r(-a  -  i>r(kHV 


Of — 1 — 3) 


as 


k  ♦  +  ®  and  a.  +  y  being  non-integer.  (If  a  +  y  is  an  Integer, 


then 


2  3 

$  (r  )  only  contains  finitely  many  terms),  and  because  a  +  ■=•  >  1  for 
n  *  ci  t 

a  >  -  y.  Thus 


11m  $  (r2)  -  1+  I  (-l)k  (*  2]  r2  -  (l-r‘) 

n-w  n,a  k-1  k 


2  *  +  7 


We  shall  estimate  now  the  difference 


9  9  ®  +  *) 

R  -  1$  (r  )  -  (1-r  )  I 

n  n  ,a 


Since 
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TT 

j-l  n+j+  j 


T(n  -  j)r(n-crt-k) 
r(n-a)f(n+k+  j) 


a  +  |  . 


holds  uniformly  in  k,  and  since  this  product  decreases  as  k  increases, 


|s  [TT(-^il) -i](-i)k(“ +  T)  r2k|  <  y+y 

k-1  j-l  n+ j+  £  k  1 


1  2 


where 


l  -  I  )cft  -  ll(*  ’  !]  r21 


k-1  j-l  n+j+  j 


N  .  .  ,  =>  |  a  +  -r  | 

<  [1  -  TT  fS=!5±J=i)]  I  (  2} 

J-l  n+j+  J  k-1  1  k  ' 


« +  i 

<  '• 


s  |rft  -iK 

k-N+1  j-l  n+j+  j 


a  + 


T,  tJk 


«  j  r^i  <  c2 -u,.. 


k-N+1  1  k 


a  +  2 


and  C^,  Cj  are  independent  of  n  and  N.  Choose  N  -  [ na+J ^ ] ,  then 
n  •*  «•  we  obtain 


as 


o+f 

1  -  <*> 


1  -  «  +  ^I7I>  '  °<V 

0+571  n 


v 


t 

I 


with  0  *  aTY/f-  ^  see  t*iat  y  an^  y  have  then  the  same  order  0(— ~ ), 

12  n° 


R  “  0(~) 

n  o 

n 


which  uniformly  holds  for  r  €  (0,1).  Therefore  we  can  write 


2  2  a  +  2  1 
$  (r  )  -  (1-r  )  +  O(i-) 

n.a  a 


and  equation  (3.4)  follows. 


Remark.  5.  In  the  case  £  »  0,  using  (2.21)  and  Lemma  1  one  can  easily  obtain 


an  -  (-1) 


,  n-1  i 

~  A 


*!•<».>♦„<«>  pn2-(lw(l» 

n 


where 


4>n(a)  " 


co 3  —  if  n  is  odd. 


if  n  is  even. 


We  shall  now  prove  the  inequalities  of  an  which  are  important  in  our 
further  error  analysis. 

Theorem  3.  If  £  <  -1  and  n>a  (a  >  -  y0, 


a  +  4 


(3.7) 


.  r  r  1  .  „  2,  2\ 

'n1  -  — T  T+  (1'r  >  ) 


ot  +  x  a  +  T’ 
2  2 
n  n 


with  the  equivalence  constant  depending  only  on  a  and  r 


ki+Zr -i 
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(Here  r  »  1  is  admissible.) 

Proof.  Since 


a 

n 


(-1)""1  C(a) 


n-cr 


-—77®  (r2)(l+0(-^) 

a+1/2  n,a  n 


holds  uniformy  in  r  as  n  ♦  ■»,  it  suffices  to  estimate 
Using  the  integral  representation  (see  [15],  p.  259): 


9 

n,a 


(r2). 


P(a,b;c;z) 


for  Re  c  >  Re  a  >  0,  larg(l-z)!  <  u ,  we  have  for  n  >  a 


(3.8) 


«n  <r2> 

n,a 


r(ff*~3/2 )  / 

r(a+3/2)r(n-a) 


tn~a-1(l-t)a+1/2(l-r2t)a+1/2 


dt. 


Writing  (l-r2t)a+1 ^2  ■  [(l-t)+(l-r2)t]a+1^2  and  noting  that  for  0  <  a 
+  i  <  1,  the  function  xa+^2  is  concave.  We  can  use  Jensen's  inequality 
([16],  p.  28)  for  a,b  >0,  r  >  s  >  0 


(3.9)  (ar+br)1/r  <  (a8+b8)1/s, 

to  obtain  for  0<a+y<l,  r  ■  1 ,  s  “  a  +  -j-: 

2a“1/2[(l-t)a+V2  +  (l-r2)a+1/2ta+V2] 

.  2^1/2(azor^±  ) 


(3.10) 


i 
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<  2^1/2f(l-^l-KL-f2)l1g+1/2  -  [(l-t)+(l-r2)t]a+1/2 

<  (l-t)a+V2  +  (l-r2)a+V2ta+1/2  . 

For  a  +  V2  >  1  the  function  xa+  ^2  is  convex,  using  once  more  the  Jensen's 
inequality  with  r  *  a  +V2  ,  s  “  1,  we  get 

(3.11)  (l-t)a+1/2  +  (l-r2)a+l/2ta+V2  <  [(l-t)+(l-r2)t]a+1/2 

-  2a“V2[(l-t)a+1/2  +  (l-r2)a+V2ta+1/2]. 

Hence  we  get  by  using  (3.8)  and  Lemma  1 

w  lihliV-  ,  I  t—-1d-t)a+V2{(l-t)a+1/2  +  (l-r2)a+V2ta+1/2}dt 

r(a+3/2)r(n-a) 

r(irt-3/2 )  rr(n-a)r(2a+2)  .  „  _2xa+  V2  T(n+  V2  )r(a+3_/2)  , 

“  r(a+3/2)r(n-a)  1  r<!tH-a+2)  Vl  r  '  r(n+a+2)  J 

22a+1r(l+a)  T(tH-3/2)  ,  „  _2xa+  V?  r(n+ V7  )r(n+3/2) 

“  r ( n+a+2 )  K  ’  r(n-a)r(n+a+2) 

/n 


(D(a) 


—  +  (l-r2)a+1/2)  (140(1)) 


a+  V? 


where 


(3.12) 


D(a)  - 


22a+1r(14-q) 

/* 


and  the  term  0(— )  la  uniform  with  respect  to  r,  thus  we  have  (using  the 
n 

notation  *■) 
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1 


(3.13) 


<5 

n,a 


(r2) 


_Qt+  V' 


+  (l-r2)a+V2 


and  Inequality  (3.7)  follows  easily.  ((3.7)  is  valid  for  r  *  1  by  the 
equation  (3.2).) 


Remark  6.  The  estimate  (3.13)  is  true  under  the  condition  n  >  a.  We  now 
consider  the  case  when  n  <  a.  In  the  proof  of  Theorem  1  we  have 


.1 .2("n)-1e 

0  (l+2r  cos  ©+r2)n~a 


d0 


2rrt"1  n! 


$ 


(2n+l ) ! !  n ,a 


(r2). 


For  n  <  a,  it  follows  that 


n,a 


(r 


2 


(2n+l) 1 1  / 
2(2n) ! !  J 


s1n2(n+1)~1e 

(l+2r  cos  9+r2  )n  a 


de 


> 


TC 

(2n+l )  ! !  ,2 
2(2n)  ! ! 


sin 


2u+l 


ode 


1 

2  * 


I 


On  the  other  hand,  apparently  we  have 


*  ,_2^  „  (2n+l)H  ,a-n  ,  2(itfl)-l„,„  .  ,a 

\,a(r  5  <  JUKTU  f  *  8ln  ede  <  4  . 


Therefore  we  have 


(3.14) 


1  $  (r2)  <  4“  if  n  <  a 

2  n.a 


and 


I  a  l  «■ 


_n-a 


a+ V2 

PI  *• 


(3.15) 


if  n  <  a 
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(see  [11],  p.  59),  we  have 
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p(a+l ,-a-l ) (cos  e) 
n 


(-1  )n  p(~“-1  »a+1)(COs(it-e)) 

a 


'(1t_e)a+1/2o(n-1/2), 


x  -  cn  ^  >  9  >  x/2, 


(3.20) 


< 

|^0(n-a-1) , 


n  >  9  >  x  -  cn 


-1 


Noticing  that  for  0  <  0  <  x,  1  -  F  =  1  -  cos  9  ■  0(9  ),  (2.15)  gives 

.a+lT 


(E> 


T ( a+n+2 ) 
O(02a+2n-a)p(a+1,-fl[-l) 


{  9a+1/2o(n"a'‘V2), 


02a+2O(n), 


cn  ^  <  9  <  x/2 


0  <  9  <  cn 


x  -  cn  >  0  >  x/2 


x  >  9  >  it  -  cn 


and  (3.16)  follows. 


4.  ERROR  ANALYSIS  OF  THE  BEST  L^APPROXIMATION  OF  (x- a)a  ON  [a,b] 

We  know  that  the  best  L2~approximation  is  given  by  the  partial  sum  of 
Legendre  expansion.  Let  Ep  ( [ a , b] ,d)  denote  the  error  of  the  best  L2~ 
approximation  on  [a,b]  of  the  function  (x-d)a,  (d  <  b)  by  polynomials 
of  degree  p.  If  [a,b]  ■  [-1,1],  then  we  will  write  simply  Ep(E)  e 
«p(  [-1,11,0. 

First,  we  consider  the  interval  [-1,1]  and  the  function  is 

u^(x)  -  (x-£)a. 
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Let 

(4.1) 

then 

(4.2) 

In  the 
-  (*-d)“ 

Lemma  2. 

where 

Now  we 
and  Ep([a 

Theorem  5. 

(4.3) 

where 


(x-E).  ~ 


I  a  p  (x), 
n=0  n  n 


on  [-1,1], 


Ep(E) 


II  y  a  P  (x)»  /  .  .  \ 

n-^fl  nn  L2  1,1) 


f  y  a2  Jl M 

‘  4  an  2n+1J 


n=»p+l 


general  case  of  the  interval  [a,b]  and  the  function  uj(x) 
the  following  relation  can  be  easily  obtained: 


Ep([a,bJ,d)  «  ha+1/2Ep(E) 


£ 


d-c 
h  * 


c 


a-tb 
2  ’ 


h 


b-a 

2  ' 


are  going  to  obtain  asymptotic  formulae  for  the  error 
,b],d). 


Ep(0 


If  £  *  -1,  then 


20c+l/2  i 

E  (~ 1 )  -  C  (a)  -  <1-K>(±» 

P  0  (p+l)2a+1  P 


(p  -*•  ») 


C0(a) 


r(l+g)  Isln  nal 


/2a+l  ii 


Consequently, 
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(4.4) 


Ep(U,b]  ,  a  ) 


c0(«) 


(p+1)2  1  P 


or 


(4.5) 


E  ([a,b] ,  a  ) 


(b-a) 

(p+1) 


a+  V2 

2a+l 


with  the  equivalence  constant  depending  on  a  but  not  p. 


Proof . 


Ep(-D 


-  {  I 

n=p+l 


j2  _L_l  v2 

n  2n+lr 


2g+1r(l+g)Z|  sin  -real  ,  r  lWd/n),  V2 

t r  1  Z  Arr4-^  J 


n=p+l  n 


4<r+3 


2g+1r(l+g)2 Isin  -real  _ 1_ 


/4a+2  (p+1) 


W  <>«<£» 


rH*.)Zl8ln  ml  2^ 1/2  . 

ir/2a+l  (p+l)Zg+1  p 


b-a 


Noting  that  h  =  (4.4)  follows  readily  by  using  Lemma  2. 


Theorem  6.  If  P  <  -1 ,  then 


(4.6) 


,(5)  -  c1(«)(^f-) 


2  a  p+1  , 

'  -rr  (1+0  (-!-)) 

(p+1)  P 


with 


a  >  0,  r 


IfJ+A2-i 


C^a) 


r(l+a)  I  sin  na  I 

/n 
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and  the  term  0(— )  Is  uniform  with  respect  to  f  <  -1 
o 
P 

Consequently, 


<*.7,  <  <!..«.«  - 


(P+1) 


with  r 


/b-d  -  / a-d 


b-a 


for  d  <  a. 


/b-d  +  /a-d  bfa-2d+2/( b-d) ( a-d) 
First  we  need  to  prove  an  auxiliary  lemma 


Lemma  3 .  If  0  <  a  <  1,  a  >  0,  then 


(A. 8) 


I  2-  “ 

_._vt  «t C7  *  s  N 

n*N  n  N 


tl  nN 


The  term  0(-^— )  is  uniform  in  s  £  10,1-e],  e  >  ( 


Proof.  Observe  that 


(4.9) 


N  ,  <*  n  n  .  . 

si  vs  Vfi  1  ^  n 

31=^  "_L  “J 


o  < 


n=N  n 


s_  1 

,,CT  1-S 


n=»N+l  Ng  nP 

N  i 
s  1 

1-s 


a-s)  7  a  -  <i£ir>a]sk 


k-l 


<  (i-s>r  1  -  fcf-)a]  7  Sk  +  (1-S)  y 

N+m  k-l  k-m+1 


[1  -  (^^(s-s^1)  +  s-1 


=  R. 


N,m" 


-  e  (c  >  0) 


(hoc—: )) 

a 


K 

s 
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Choose  m  =  ,  then  for  N  ■*  *  we  obtain 


*n7 


V-  4  1  -  (1  +  ?>'0  +  s” 


-  —  P  S  +  i-  +  q((^—  — 

M  „  1  N  UU  N  ;  ; 

N  In  — 
s 


0(^~ .) 


and  this  term  0(— - — )  holds  uniformly  with  respect  to  0<s<1-e, 
e  >  0.  1 


Proof  of  Theorem  6.  Lemma  1  and  Lemma  2  give 


V5> 


{  I 

n=p+l 


,2  _J — iV2 

n  2n+l ' 


"(1+g)  1  sin  Ttgl  (l-r2)a+ 


/it 


(2  r)c 


{  l 

n=p+l  n 


-2n  „  V2 


2a+l  2n+l 


}  (1  +  0(~)> 


..  2.a+Vo  ,  2(p+l)  l/o  , 

(a)  ILzlJ - {*  _E — —A  (1  +  0(^£))(1  +  0(~)) 

1  (2r)a  1-r2  ( p+1 )2a+2 ’  P  p° 


1-r2  a  r1*1  1 

ci(a)(lT)  7~T£+r  (l  +  °(ig)) 

(p+i)  p 


a+  V2 


1 


where  a  *  -  €  (0,1)  and  the  term  0( — )  is  uniform  with  respect 

g+3/2  pCT 

to  0  <  r  <  1  -  e,  e  >  0.  Considering  the  general  case  we  obtain 


d-c  _  2d-a-b 
h  “  b-a  * 


b-a 


/b-d  -  /a-d 


1 1 1  +Vt^l  bfa-2d+2/(  tr-d)  (a-d)  /t-d  +  /a-d 
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(a  <  d) ,  and  (4.7)  easily  follows.  " 

Remark  7 .  If  d  <  a,  the  singular  point  Is  outside  of  the  Interval 
[a,bl,  then  the  error  Ep([a,bl,d)  reduces  exponentially  by  the  above 
theorem,  and  the  rate  Is  characterized  by  the  ratio  r.  This  value 
depends  only  on  the  ratio  of  the  length  of  the  interval  and  the  distance 
between  the  position  of  the  singularity  and  the  Interval,  but  is 
independent  of  a.  In  Section  5  we  will  see  that  this  property  also  holds 
for  more  general  functions.  In  fact,  let  III  =  b  -  a,  6  “  a  -  d  = 
dist.  (d,[a,b]),  then 


=  -  (1  +  lUid).) 

n  D-a 


-  (1+2X) 


where  X  **  fyp  Therefore, 


(4.10) 


1+2X+2  A(l+X) 


Thus  r  is  geometric  invariant  of  the  error. 

We  are  now  going  to  obtain  the  estimate  which  is  uniform  with  respect 
to  E  <  -1.  We  need  first  the  following  lemma: 

Lemma  4 .  Let  0<s  <1,  o>0,  N  >  2. 

1)  If  0  <  s  <  1  "  i,  then 


(4.11) 


n“N  n 


®  n  N  , 

y  §_  .  s__i 

_<J  „<J  l-s 


2)  if  1  -  <  s  <  1, 


then 
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N  N 


(4.12) 


„  o+l  „o 

n"N  n  N 


These  Inequalities  hold  uniformly  in  s  and  N. 


Proof.  As  in  the  proof  of  Lemma  3,  we  have 


N  ,  »  n 

s  1  y  8 

s  1 
„o  1-s 


If  0<s<l<~— ,  then  for  any  m  >  1 


Choose  m  »  N  and  notice  that  (1  -  jj-)N  ♦  we  get 


thus 


N  ,  “  n 

S  1  _  r 

Mo  1-s  „  o  ,  , 

0  <  5 - - — --N-—  <  1  -  —  (1  -  -), 

s_  _1_  2°  e 

„o  1-s 


N  i 
8  1 


oo  n 

>  i  a- 

L  o 


N  1-s  n”N  n 


,  ,  N  , 

>  i-  (1  -  I)  2-  -1_ 

-o  e  a  1-s 


and  (4.11)  follows. 

In  order  to  show  (4.12),  observe  that 


T  <  aN  T  1  t  .N  C  dt 

Li  0+1  Li  0+1  1  1  0+1 

n^I  n  n«N  n  N-l  t 
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If  1  >  s  > 


since  (1  - 

This  proves 

Theorem  7. 
1)  if 


s  1 


0  (N  — 1 )  ' 


N  »  i 

S  1_ 

a  «~1 '  Na 


2°  sN 

<  —  —  for  N  )  2. 

a  N0 


1  -  — ,  then 
N 


®  X 


oo  n  * 

>  -fe  >  ^  dx 

*N  n  N  x 


I 


0  (y+N) 


a+1 


dy 


(y  -  x-N ) 


si  r  a-l/N)y  dv 

N°  0  (1+y/N)0  N+y 


N  »  ,,  .  /H.Nt 

2-f 

N0  0  (l+t)^1 


(t  -  y/N) 


1/N)N  t  i  as  N  -  «,  (1  -  I)N  >  ^  for  N  >  2.  Thus 


n«N  n 


o+l  .,o  J 


sN  4-t 


IT  0  (l+t) 


o+l 


dt. 


(4.12). 


Let  t  <  -1,  p  +  1  >  a,  then 


2  1 

0  <  r  <  1  -  —rr,  then 

p+1 


p+l-a 


vs)  -  ■  ,a+u 


r- 


1  +  (i-r2)a+V2] 


/1-r2  (pfl)  (p+l)a+  '2 


(4.13) 
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2)  tf  1 - 77-  <  <  1  then 

p+1 


(4.14) 


E  (E)  -  -r  — -  3-  r - 1 - -  +  (1-r2)a+1/2l 

(p+l)a+  72  (p+l)a+  <2 


Inequalities  (4.13)  and  (4.14)  hold  uniformly  in  r  and  p,  where 


1 


IeI+/e 


In  the  case  of  general  interval  [a,b]  the  inequalities  have  to  be 
modified  by  multiplication  with  a  factor  (b-a)®*^^ 


Proof.  The  results  are  the  direct  consequence  of  Theorem  3,  Lemma  4,  and 
the  simple  inequality 


•77  (/x  +  /y)  c  /x+y  <  /x  +  /y.  ■ 

/2 


Remark  8 .  By  Remark  6  these  estimates  are  also  valid  for  pi-1  <  a. 

Now  we  shall  explore  the  error  behavior  for  the  case  -1  <  E  <  1.  As 
we  see  in  Theorem  2  the  coefficients  an  behave  in  a  more  complicated 
manner,  hence  we  cannot  obtain  a  simple  asymptotic  formula  as  in  the  other 
cases.  However,  we  have  the  following  estimates: 


Theorem  8.  For  -1  <  E  <  1,  there  exist  a  constant  C  >  0,  which 
depends  only  on  a  such  that 


(4.15) 


Ep(0  < 


0  <  9  <  x  - 


p+1* 


n  -  —  c  9  <  x, 
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where  6  ■  min(6,w-0),  9  »  arccos  £. 

For  the  error  in  the  general  interval  [a,b],  a  <  d  <  b,  the  right- 
hand  side  has  to  be  multiplied  by  (b-a)orh^2>  and 


0  ■  arccos 


2d-a-b 

b-a 


Proof.  By  Theorem  4  we  can  write 


C  - 
n 


2o+T 


*n  I  <  \ 


7T  —  —  <  8  <  TT 

n 


0  <  9  <  TT - 

n 


with  C  -  C(a),  6  •  min(TT-0,0). 

Therefore,  if  ir  -  <  0  <  tt,  then 


«.«»  -  ( i 


n»p+l 


2p 

‘  Cl  !•  -4^3  +  1 


l/2 

_E _ i  z 


n**p+l  n  n“2p+l  n 


2a+2 


<  C(^r)20*1. 


The  case  when  0  <  9  <  r  -  follows  similarly. 

For  the  other  side  of  the  inequality  we  have  the  following  result. 

Theorem  9.  Let  e  <  9  <  n  -  e,  e  >  0,  then  there  is  C  ■  (a,e)  >  0 
such  that 


(4.16) 


V"  ’  V 


6  vor+-  V2 


(6  "  min(9,n-0)  ) 


the  estimate  is  uniform  with  respect  to  6  (5  ■  arccos  £)• 
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T 


Proof-  By  Theorem  1 


(4.17) 


EP(E)  '  /|r(l+a)(sin  0)a+1/2{  l  J^  (a  +  T>^]  +  oC-^)}1^ 

n-p+1  aia+z  n2a+3  ’ 


*  »-U  »2«+2  1  f7^h 


The  term  0(  a+^)  Is  uniform  with  respect  to  9  E  [e,n-e]. 

p 

It  can  be  readily  seen  that  (4.16)  Is  proved  if  we  show  that 


(4.18) 


max  {  cos2  [ (m  +  1)0  -  (a  +  |)^]}  >  sin2  4. 


m-n.n+l 


We  prove  this  by  contradiction.  Let  y  -  -|,  thus  sin  |  - 
cos(y  -  y  y).  If  (4.18)  is  not  true,  then  for  m  »  n  and  m  »  n  +  1  we 


cos2[(m  +  |)0  -  (a  +  |)£]  <  Cos2(|-Yy). 


Thus  there  are  Integers  k(n)  and  k(n+l)  such  that 


Therefore 


-Y  <  (2n+l)|  -  (a  +  i+  2k(n))  <  Y , 
-y  <  (2n+3)|- -  (a  +  y  +  2k(n+l))  <  Y  > 


y  <  J  -  (k(n+l )  -  k(n) )  <  y . 


I  * 


I  *  •* 


0  9 

because  y  *  —  and  y  <  1 - ,  we  get 

it  it 


k( n+1 )  -  k( n)  >  —  -  y  >  —  -  —  *  0, 

7i  it  n 


k(iH-l)  -  k(n)  <  —  +  y  <  i+i-I  -  1. 

it  it  it 


Since  k(n+l)  -  k(n)  Is  an  Integer,  we  have  the  desired  contradiction. 


Remark  9.  It  Is  easy  to  obtain  the  estimate  for  the  case  £  *  0: 


(A. 19) 


Eq(0)  a 


C(a) 

(pfl)a+  V2 


where 


In  fact  by  (3.6)  we  have 


r(n-g) 

/(2a+l  )it 


Ep(°) 


A  (  L,  *”(',>2 

n-p+1  n 


/-  »  (sin  -|2-)2  =>  (.cos  -T- > 

\  rd+a)  {  i  — - hr+  ?  — 

^[£+1]  (2k)2o(+2  ,  (2k+l)2a+2 


(cos  V2 


5  /!r(14*>  hk  ? 


_i _ .1V2 


2  k-[* k 


s  rd-hz)  1 

/(2a+l)tr  ({y4.1)a+1/2 


At  last  we  consider  the  problem  of  the  asymptotic  behavior  of 
Ep([a,b],^)  as  the  size  III  of  the  Interval  approaches  zero.  For 
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I  '  ■’ 


4V.)  _  JL|  <  3 1 1 1  3  ,  , 

m  4x|  4e(4e-3E)  “  J 

4e 


|l|™0  I1 t  "  4x  unlformly  hold  with  respect  to  x  >  e,  E  >  o. 


By  (2.16) 


a  =  C_13n~l  r(l+a)sln  Tig  r(n-<z)  n-a  .  ? 

n  _  ~  ‘ —  r  $  (r  ). 

/n  •  2a  r(n+l/2)  n,a 


Recall  that  $n^(r2)  has  a  majorant  series  J  — „  ( 


k=l  k 


a+ 3/2  <  00  lsee  proof 


of  Theorem  2),  thus  11m  $  (r  )  »  1  i  , 

„  n,au  ;  1  holds  uniformly  in  n.  Therefore 


we  have 


„,L  *«  k?t  <  ««>  ?  r-itol  ]2  2  r2'»-«> 

n  P+2  n*p+2  r(n+V2)  2n  1 


<  c(a)  ~~  r2(p+2-a)( 

1-r 


it  follows  that  If  X  >  e  and  |l|  is  small,  then 


i2  *  rili^2a:+l  ft  22 


Vl(I’0)  *  {  I 


n=p+2 


n  2n+l ' 


||2a+12(p+2-0[) 
<  C(a)  — - L_ - 


<  C(a,c) ll |2p+  . 


Hence 


2 

f  — 2— - — .) 

Ill  ^3/2 


_1 _  <-  ll  I  -i2ct+l  ,  2  2  ? 

,2p+3  2  1  *  fVl  2PT+  Vl(1>0) 


i 
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Theorem  11.  Let  u  €  l^(a,b) 


satisfy  the  above  condition  (5.1),  and  let 


I 

h 


[a,b], 


d  =  dist.(xQ,I) 

b-a 
2  ' 


Then  the  error  of  best  I^-approximation  by  polynomials  of  degree  p  can 
be  estimated  as  follows:  if  e  <  r  <  1  -  e  (e  >  0),  then 

.  2  a 

(5.2)  E  [a.b]  <  KCCcOh0^  '2  r1*1 

P  2  r 

where 

r  5  r(^)  “  - 1 

1  +  i  +  /f  <2  + 

h  h  h 

Proof.  First  of  all,  we  make  the  linear  transformation  as  before: 


x  =  c  +  ht,  t  f  [-1,1] 


c 


a+b 
2  ’ 


h  - 


b-a 

2  • 


It  maps  [-1,1]  onto  [a,b].  Therefore  for  the  function 


w(t)  =  u(c+ht) 
defined  on  [-1,1],  the  singular  point  will  be 


6  ?  dist . (Eq  ,  [-1 , 1  ] ) 


1 


and 


I  »  ' 


=  dist.(x0,  [a  ,b] )  - 


and  the  growth  condition  (5.1)  becomes 


Iw(c)  -  u0l  t  K|(c+hc)  -  (c+hC0)l“ 


K  hak  -  50la. 


As  usual,  we  expand  w(t)  int  the  Legendre  series 


n=0 


(t  e  [-i,i]) 


(5.4) 


Ep[-l,U 


1  1  a2  — _ fa 

n-U  n  2n+1 


Without  loss  of  generality,  we  may  assume  Xq  <  a,  thus  £g  <  -1  and 

?g  =  -1  ~  6. 

Since  w(t)  is  analytic  except  at  f  [-1,1],  we  have 


an  “  /  w(t)Pn(t)dt 


2|±i_L_J  w(t)[t2-l)n](n)  dt 
2  2n-n!  -1 


dt 

2  2n«  n!  -1 


2-f --±-1  (fk*  -^itla-t2)"dt 

2  2n-n!  -1  21,1  Y  (C-t)"*1 
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1 


where  y  is  any  contour  to  which  t  is  interior  and  rg  is  exterior, 
and  it  is  positively  directed. 

Now  we  choose  y  to  be  the  circle  centered  at  t  with  a  radius  R  * 
R(t)  <  It  -  ?^q  I  ,  then  we  have  the  following  estimates; 


w(C) 


2711  Y  (C-t)n+1 


.  w(C)-ur. 

dcl  =  *27ti  ^  ,,  TTn+r  dc' 


Y  (C~t ) 
lw(C)-u 


0 


.  rrt-1 


ds 


2u 


<  »L  (-1  f  |C  -  C 

n  l2it  >  ^0 


d9l 


where  C  =  t  +  Re*'®. 


Let  be  the  intersection  of  the  circle  to  the  segment  [Eg*1!* 

and  let  o  »  f  -  be  the  distance  of  P.q  from  the  circle,  then 

2  it 

-  —  /  1C  -  ?nla  d9 

2  it  0 

1  rn  r  2  9 . 2  ,  »  _  2  0. a/2 

m  Jn  J  fP  +  ^2R  cos  j)  +  4PR  cos  J 1  d®- 

-i n 


i 


«  i 
a 


M(o) 


I  I  • 


If  -  —  <  a  <  0,  the  Integral  Is  dominated  by 


I  2R  cos  j|  a  , 


and  if  a  >  0,  it  is  continuous  in  o.  If  we  set  RQ”t-fQ”R+o, 
we  can  let  p  -*•  0+  so  that  R  *  Rq,  and  obtain 


M  =  lim  M(p ) 
p*(T 


[  !2R  cos  |l“  d9  -  C(ot)Rp 


Because 


C(«)  -  ~  I  \2  cos  ||a  d9. 

-it 


.  it  ,oc+l  u/2 

—  J  |2  cos  dQ  *  -  J  (sin  $)“  dij), 


using  the  fact  that  for  0  <  $  <  — , 


we  obtain 


—  <J)  <  sin  i>  <  i>, 

'll 


0a+l  -a  ,  c H 

2  2  1  n 

tt  tt  a+1  2 


<  C(a )  < 


,a+l  ,  a+1  _a 

l  l  71  Tt 

u  a+1  2  a+1 


this  holds  for  all  a  >  -1 .  Therefore  we  have 


la  I  <  2£±L./‘H^<1-tVdt 

”  2  2“  -l 
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KC(a)ha  ^±1  /  _g-t2)n,  d 

2n+1  -1  ( t+1+6 )n_a 


In  Section  2  we  obtained 


/  -^I-t2  dt 
-1  (t+l+6)n  a 


(2  r')n-a  2trhl  n!  *  /  2, 

( 2 n+1 ) ! !  *n,afr  } 


with  $n  ^  defined  by  (2.8),  and 


l+6+/(l+6)2-l  l+6V6(2+6) 


thus  (5.5)  gives 


‘  KC(«,h“(2'')”"“  (JJTJTTT  »„ia<r2> 


Now  all  estimates  of  ^  obtained  in  Section  3  can  be  applied  here, 
for  example,  we  have 


®n,«"2)  ‘  <l-r2>*+V2  *  O(i-) 

n 


with  0  <  r  <  1,  a  >  0,  uniformly  hold  in  r.  Since  we  also  have  the 


asymptotic  equality 


it  follows  that 


(2n)l!  _  /to 

(2n-l ) !  !  =  v  2 


lanl  <  KC(a)hs  /n  r"'“(l-r2)I+l/2  , 


At  last,  we  obtain  for  e  <  r  <  1  -  e  (e>0): 


r: 

► 


Epf-l,1]' 


2  2 


n»p+l 


n  2n+l 


and 


<  Y  [KC(a)]2  h2a(l-r2)23+1r2(n'a) 


[KCCa^riZlil0  r^1]2 


SpU.b]  -  /2hE  1-1,1] 


<  KCCalh^fiZlV  r^1. 


Remark  11.  As  In  the  case  discussed  in  previous  sections,  If  Xg  is  an 
endpoint,  then 

ha+V2 


E  [a,b]  <  C 


(p+1) 


2a+l 


and  if  Xg  Is  interior  to  (a,b],  then 


a+  V2 


E  [a,b]  <  CMU 

P  P+1 


Finally,  the  case  when  Xg  is  not  real  also  gives  an  exponential  rate  of 
convergence: 

E  [a,b]  <  Cha+  ^2  rp. 

P 

Here  r  is  determined  as  follows:  let  pg  be  the  sum  of  semi-axes  :>f 

the  ellipse  T)  which  has  the  foci  at  ±  1  and  passes  through  the  point 
P° 

rQ  -  with  c  -  h  -  then 


r 

I* 
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po  "  T  { lso  _l1  +  uo  +  11 

and  —  <  r  <  1.  This  result  follows 
p0 

p.  245. 


+  /<Ie0  -ll  +  Ie0  +  il>2  -  4  | 

easily  from  Theorem  9.1.1,  of  [11), 


6.  ERROR  ANALYSIS  OF  p-VERSION  OF  FEM  FOR  THE  MODEL  PROBLEM 

Since  the  error  of  the  finite  element  solution  u  of  the  model 
problem  is  exactly  the  error  of  the  best  Lj-approximation  of  u'  by 
piecewise  polynomials  (see  Section  1),  the  results  of  previous  sections 
give  the  error  analysis  for  the  p  version  of  FEM  with  only  one 

Cl“l 

element.  One  only  needs  to  notice  that  u'(x)  *  a(x~E)+  -  const,  so 

that  the  error  estimates  will  be  obtained  by  replacing  a  in  the  previous 
results  by  a-1,  and  p  by  p-l,  and  taking  into  account  the  change  of 
length  of  the  interval.  Thus  the  following  results  follow  easily  from 
Theorem  5,  6,  7,  8  and  9. 

Theorem  12.  Let  Ep(E )  be  the  error  of  the  finite  element  solution  of 
the  model  problem  (1.1)  when  using  only  one  element  I  *  [0,1]  itself. 

(E  is  the  position  of  the  singularity.)  Then  for  a  >  y,  one  has 

1)  if  E  -  0,  then 


V°>  ’  °0<a>  -2arT  +  0(?» 


(P  +  00  ) 


with  Cg(a) 


ar (a)  lain  »g  I 
x/ 2a-l 


2)  if  E  <  0,  then 


V°  '  ci<0)(-4r)°'1  4ri +0't» 


mesh.  As  p  +  °°  the  rate  of  convergence  will  be  the  same  as  in  the 

case  £  «  0,  i.e.,  — 2~[ *  5  is  in  the  interior  of  any  mesh 

P 

interval,  then  the  rate  of  convergence  for  p  <=  will  be  the  same  as 
when  0  <  £  <  1  and  it  is  half  in  the  exponent  of  the  case  above,  i.e. 


large.  For  samll  p,  the  relation  of  the  error  in  dependence  on  the 
number  of  total  degrees  of  freedom  of  the  finite  element  space  has  in 
general  two  phases.  If  the  mesh  is  properly  designed  (in  general,  there 
is  sufficiently  strong  refinment  around  the  singularity),  then  the  rate  of 
error  reduction  will  be  exponential  in  the  beginning  phase.  When  p  is 
large  enough,  this  rate  becomes  algebraic.  If  the  refinement  of  the  mesh 
is  not  strong  enough  (for  example,  the  uniform  mesh),  then  the  exponential 
part  of  the  error  reduction  cannot  even  appear.  If  £  lies  outside  of 
the  interval  then  (also  for  one  element  mesh)  only  the  exponential  phase 
appears.  Typically,  when  £  lies  inside  the  interval  the  graph  of  the 
error  in  dependence  on  N  in  double  logarithmic  scale  is  S-shaped 
(actually  reflected  S).  The  first  part  is  the  exponential  phase  and  the 
second  the  algebraic  one  (see  Fig.  7.2).  In  practice  we  would  like  to 
achieve  desired  accuracy  in  the  exponential  phase.  Using  Section  5  the 
results  can  be  extended  to  the  general  case  of  a  function  with  a 
singularity  of  type  xa. 

7.  NUMERICAL  RESULTS 

In  the  previous  section  we  have  shown  various  estimates  characterizing 
the  error  behavior  of  the  p-version  of  FEM.  Here  we  will  numerically 
analyze  the  accuracy  of  the  estimates,  the  range  of  asymptotic  validity, 
etc.  The  error  in  this  chapter  is  measured  in  the  energy  norm.  As  said 
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in  Section  6,  the  estimates  are  obtained  from  the  estimates  of  the  error 

in  l^-norm  by  replacing  a  by  a-1  and  p  by  p-1 . 

In  the  tables,  p  is  the  polynomial  degree  of  the  approximation,  and 

the  error  of  the  finite  element  solution  (when  only  one  element  is  used) 

A  A 

is  denoted  by  E„  =  E  (a.C).  Let  E  -  E  (a.O  be  the  error  given  by 
P  p  p  P 

B  B 

certain  asymptotic  formulae,  E  •  E  (ot.O  be  given  by  some  some 

P  P 

estimates-  We  will  compute  the  ratios 


R 


A 

P 


These  ratios  will  be  called  numerical  constants  which  reflect  the  quality 
of  the  estimates,  the  range  of  the  asymptotic  validity,  etc. 

In  the  case  C  “  0,  the  asymptotic  formula  is  (cf.  Theorem  12) 


(7.1) 

where 


and  we  have 


CQ(a) 

2«-l 

P 


CQ(a) 


aT(a)  isln  »a I 
ir/2a-l 


lim  R 

p+<» 


A 

P 


1. 


The  numerical  results  are  shown  by  Table  7.1  for  a  *  0.7  and  a  ■ 
3.5.  We  can  see  that  for  a  small  a  (the  singulairty  is  strong)  formula 
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TABLE  7.1 
£  -  0. 


p 

a  * 

0.7 

a  *  3. 

5 

E 

Rn 

E 

Ra 

P 

P 

P 

P 

1 

4.743  E-l 

0.9877 

1.021 

0.2032 

2 

3.627  E-l 

0.9967 

3.402  E-l 

4.335 

3 

3.090  E-l 

0.9985 

3.093  E-2 

4.488 

4 

2.756  E-l 

0.9992 

2.379  E-3 

1.940 

5 

2.522  E-l 

0.9995 

4.760  E-4 

1.480 

6 

2.344  E-l 

0.9996 

1.400  E-4 

1.300 

7 

2.204  E-l 

0.9997 

5.154  E-5 

1.208 

8 

2.090  E-l 

0.9998 

2.210  E-5 

1.153 

9 

1.994  E-l 

0.9998 

1.057  E-5 

1.118 

10 

1.912  E-l 

0.9999 

5.495  E-6 

1.094 

11 

1.840  E-l 

0.9999 

3.053  E-6 

1.077 

12 

1.777  E-l 

0.9999 

1.790  E-6 

1.064 

13 

1.722  E-l 

1.000 

8.999  E-7 

1.054 

14 

1.671  E-l 

1.000 

6.978  E-7 

1.046 

15 

1.626  E-l 

1.000 

_ 

4.585  E-7 

1.040 

(7.1)  gives  very  good  results  even  for  small  p.  Thus  In  this  case  the 
asymptotic  range  Is  quite  large.  For  a  large  (the  singularity  Is  weak 
or  the  function  Is  "smoother"),  the  asymptotic  range  shift  to  large  p, 
but  for  small  p  the  accuracy  of  the  asymptotic  formula  Is  still  quite 
good.  It  Is  also  seen  that  there  is  a  big  gain  of  the  error  reduction 
when  p  Increases  from  1  to  2,  3,  especially  for  large  a. 

In  the  case  £  <  0,  the  error  reduces  exponentially.  By  Theorem  12 
we  have  the  following  asymptotic  formula: 


(7.2) 


.  2  a-1 

-  V-X^) 


p 


C^a) 


ar(a)|gln  wa| 


where 


I  >  •' 


/1-r  -  /-y 

/I^F  + 


lim  R  =  1. 

P~  P 


In  this  case  we  also  have  the  estimates 


C.(a)  ptl-a  .  _  a-V2 

J -  I -  (1-r2)  ] 

£7  ,«-■* 


2  1 
If  0  <  r  <  1  -  - 
P 


C2(a) 


p+l-a  Df  v  2  a-V2 

- -  f£i£2_+  (1-0 


a-  V2  a~  V2 

P  L  P  1 


If  1  -  -  <  r2  <  1 
P 


C2(a)  -  C1(a)/2a_1 


22g~1r(a) 


Here,  D(a)  Is  the  number  obtained  In  proving  Theorem  3.  (See  (3.12)* 

C2(o)  Is  chosen  so  that  as  p  ■*  ®,  the  first  formula  asymptotically 

agrees  with  (7.2),  therefore  we  also  have  lim  R  -  1.  It  is  shown  in 

n  IITI  P 
P*® 

Table  7.2(a)-7.2(c)  that  the  asymptotic  behavior  Is  not  too  simple  (it 

will  be  described  later)*  Note  the  two  parts  of  the  formula  E  coincide 

P 


i  » 
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at 


2 


r  « 


differ. 


1  -  — .  For  r  “  1  (l.e.  £  ”  0)  the  formulae  (7.2)  and  (7.3) 

P 

In  fact,  if  we  write  (7.2)  by  EA(0),  then 

P 


EB(0)  _ 

-2 -  =  J  4a-2. 

EA(0) 

P 


The  ratio  R^  in  this  case  is  bounded  above  and  below  by  constants  which 
depends  on  a. 

It  can  be  seen  in  Table  7.2  that  the  asymptotic  range  of  (7.2)  depends 

on  a  and  £.  When  a  increases  or  £  gets  close  to  the  approximation 

interval  [0,1],  the  asymptotics  are  shifter  toward  large  p.  If  a  is 

A  B 

large  and  £  is  close  to  -1,  E^  is  large,  while  E^  gives  good 
estimation  in  the  sense  that  for  large  range  of  both  o  and  £  the  ratio 


Rp  is  quite  stable. 


TABLE  7.2(a) 

£  -  -0.0005  (r  -  0.9563) 


p 

a  -  0.7 

a  -  3.5 

L 

E 

RA 

rb 

E 

ra 

rb 

P 

P 

P 

P 

P 

P 

1 

4.192  E-l 

0.4780 

0.6334 

1.022 

3078. 

0.0549 

2 

2.945  E-l 

0.5705 

0.5806 

3.404  E-l 

12130. 

1.224 

3 

2.323  E-l 

0.6251 

0.5440 

3.091  E-2 

4761. 

1.324 

4 

1.928  E-l 

0.6636 

0.5162 

2.373  E-3 

1047. 

0.5972 

5 

1.647  E-l 

0.6929 

0.4939 

4.738  E-4 

476.9 

0.4753 

6 

1.433  E-l 

0.7165 

0.4752 

1.390  E-4 

276.9 

0.4352 

7 

1.268  E-l 

0.7359 

0.4593 

5.106  E-5 

182.4 

0.4212 

8 

1.126  E-l 

0.7524 

0.4453 

2.180  E-5 

130.0 

0.4187 

9 

1.010  E-l 

0.7666 

0.4329 

1.038  E-5 

97.80 

0.4224 

10 

9.116  E-2 

0.7791 

0.4218 

5.373  E-6 

76.56 

0.4297 

a  =  0.7 


TABLE  7.2(b) 

0.05  (r  -  0.6417) 


2.190  E-l 
9.671  E-2 
5.063  E-2 
2.640  E-2 
1.475  E-2 
8.434  E-3 
4.904  E-3 
2.886  E-3 
1.716  E-3 
1.027  E-3 


0.7478 

0.8362 

0.8765 

0.9003 

0.9162 

0.9276 

0.9362 

0.9429 

0.9484 

0.9528 


0.4700 

0.4321 

0.4706 

0.4962 

0.5150 

0.5298 

0.5417 

0.5518 

0.5604 

0.5678 


1.153  | 

3.580 

E-l 

2.912 

E-2 

1.937 

E-3 

3.254 

E-4 

7.844 

E-5 

2.322 

E-5 

7.853 

E-6 

2.920 

E-6 

1.165 

E-6 

0.0340 
129 
934 
227 


3.729  E-2 
4.302  E-3 
5.744  E-4 


TABLE  7.2(c) 

-1  (r  -  0.1716) 


a  -  0.7 

a 

-  3.5 

E 

ra 

rb 

E 

ra 

Rn 

P 

P 

P 

P 

— 

P 

P 

33 

4.716 

65 

5.981 

15 

1.696 

0.0191 
1.19 
2.42 
1.58 
1.54 
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For  the  case  0  <  £  <  l,  we  use  the  estimate: 


(7.4) 


'1-5 


2  a  -Vo 


0  <  9  <  it  -  - 

P 


,2a-l 


71  -  —  <  0  <  TT 

P 


where  9  ■  arc  cos  ?.  In  this  case  we  cannot  have  an  accurate  asymptotic 

B  B 

formula  as  for  £  <  0.  The  numerical  constant  R  *  E  /E  no  longer 

p  p  p  6 

shows  a  monotonic  behaviour.  Tables  7.3(a)(b)  are  made  for  a  *  .7  and 
3.5  and  ?  ranging  from  0.0005  to  .75.  The  numerical  results  show 

g 

that  the  ratio  R  is  stable  (it  is  bounded  above  and  below  by  constants 
P 

g 

which  depend  on  a).  For  small  a,  the  range  of  R  is  quite  small 

g 

(this  means  that  E  gives  very  good  estimation).  When  a  is  large, 

P 

this  estimate  is  not  so  good  but  we  still  have  the  right  rate  of 

g 

convergence.  An  interesting  fact  is  that  the  error  E^  for  fixed  p  is 

nearly  symmetric  in  £  (about  C  ”  ),  yet  the  function  (x-£,)+  is 

not  symmetric.  We  also  see  that  the  error  (for  fixed  p)  has  the  maximum 

at  5  *  V2  (the  middle  of  the  interval).  It  tends  to  be  smaller  when 

5  moves  to  the  endpoint  of  the  interval. 

The  graph  of  the  error  function  E  (a.O  for  a  »  1.5  is  shown  in 

P 

Figure  7.1.  It  is  very  clear  that  if  the  singularity  is  located  outside 
of  the  Interval  of  approximation,  then  the  rate  of  convergence  is 
remarkably  increased.  This  fact  is  important  in  designing  a  right  mesh 
for  the  FEM.  It  will  be  shown  in  the  second  part  of  the  paper  that  a 
strong  refinement  around  the  singularity  will  greatly  reduce  the  error  for 
the  same  number  of  degrees  of  freedom. 


1  I  •' 


TABLE  7.3(a) 
a  =  0.7. 
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TABLE  7.3(b) 
a  »  3.5 


t  -  0. 

005 

E  =  0. 

0005 

E„ 

rb 

f 

rb 

P 

P 

p 

P 

1.008 

2.304 

1.019 

3.384  E-l 

49.50 

3.400  E-l 

3.114  E-2 

51.88 

3.095  E-2 

2.436  E-3 

22.81 

2.424  E-3 

4.983  E-4 

17.80 

4.778  E-4 

17.07 

1.508  E-4 

16.08 

1.410  E-4 

15.03 

5.749  E-5 

15.46 

5.211  E-5 

14.01 

2.566  E-5 

10.70 

2.241  E-5 

13.43 

1.287  E-5 

7.640 

1.076  E-5 

13.07 

7.070  E-6 

5.758 

5.622  E-6 

12.85 

E  «  0 

.5 

E  -  0 

.75 

E„ 

R® 

E 

rb 

P 

P 

P 

P 

1.552  E-l 

0.3548 

2.091  E-2 

0.0736 

9.953  E-2 

1.820 

1.711  E-2 

0.4817 

4.289  E-2 

2.647 

1.202  E-2 

1.142 

8.816  E-3 

1.290 

6.987  E-3 

1.574 

3.630  E-3 

1.037 

3.097  E-3 

1.362 

1.897  E-3 

0.9364 

9.555  E-4 

0.7263 

1.128  E-3 

0.8843 

6.864  E-4 

0.8283 

7.293  E-4 

0.8534 

5.736  E-4 

1.033 

5.001  E-4 

0.8334 

2.932  E-4 

0.7521 

3.901  E-4 

0.8196 

2.127  E-4 

0.7485 
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tog  Ep{*) 


Figure  7.1.  p-version  with  one  element  (a  “  1.5). 

So  far  we  discussed  the  p-version  with  one  element  only  (m  "  1). 
Let  us  now  address  briefly  the  case  when  the  number  of  elements  m  >  2. 
We  assume  that  the  solution  is  u(x)  »  xa  -x  (i.e.,  the  case  K  “  0). 
Figures  7.2(a)  and  7.2(b)  show  the  cases  when  a  ”  0.7  and  a  *  l.l, 
respectively.  In  the  figures,  m  »  2  ~  10,  and  the  meshes  are  made  bv 
geometric  progression  with  a  ratio  q  ■  0.15.  The  dotted  curve  is  the 
error  for  the  nearly  optimal  mesh-degree  combination  (cf.  Part  2).  In 
Figure  7.3  we  compare  the  cases  in  p-verslon  for  uniform  mesh  and 
geometric  mesh  (a  ■  0.7). 


Figure  7.3  (a  m  0.7). 

a:  uniform  mesh  (m  *  50); 
b:  uniform  mesh  (m  *  100) 
c:  geometric  mesh  (m  ■  20,  q  -  0.15). 

It  is  clearly  seen  that  the  p-version  has  often  two  phases.  The 
first  phase  (when  p  is  small)  has  approximately  an  exponential  rate. 
The  second  phase  (when  p  is  large)  tends  to  be  algebraic.  This  figure 
will  be  called  an  S-curve  (actually,  it  is  a  reflected  S).  For  the 
geometric  meshes,  whan  having  the  same  number  of  elements,  if  q  is 
small,  the  range  of  exponential  phase  is  enlarged.  However,  it  does  not 
mean  that  the  smaller  q  is  the  better  the  result,  since  the  curve  is 
also  shifted  up.  Roughly  speaking,  when  a  is  small  or  the  required 
accuracy  is  high,  it  is  better  to  use  smaller  q.  When  a  is  large  or 
required  accuracy  is  low,  large  q  is  preferable.  But  the  difference 
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for  q  «  0.1  ~  0.3  is  not  too  large  if  the  required  accuracy  is  not  too 
high  (say  12  ~  52) . 

Secondly,  for  a  small  (solution  is  highly  unsmooth)  one  needs  more 
intervals  than  in  the  case  of  a  large  (solution  smooth).  For  example, 
when  a  ■  0.7,  using  m  m  2  or  m  ”  5  one  even  cannot  obtain  the 
accuraacy  of  52  in  practice  since  the  polynomial  degree  cannot  go  too  high 
(say,  p  ■  10).  Figures  7.2(a),  (b)  show  all  these  curves  lie  above  the 
curve  for  the  optimal  h-p  extension  (see  also  Part  2).  Likely  there  is 
an  envelope  for  these  p-version  curves,  and  this  envelope  lies  above  the 
optimal  h-p  extension  curve. 

It  is  very  clear  that  the  mesh-design  is  crucial  to  achieve  a  good 
rate  of  convergence.  If  the  singularity  is  present,  then  the  uniform  mesh 
with  few  elements  is  not  acceptable.  Figure  7.3  shows  that  the  uniform 
mesh  performs  badly.  There  is  no  "exponential"  phase  at  all  and  for  a  » 
0.7  it  cannot  get  even  52  accuracy  for  hundreds  of  times  many  degrees  of 
freedom  comparing  with  the  case  when  a  geometric  mesh  (q  ■  0.15,  m  * 

20)  was  used. 

In  practice  (see  [2])  one  would  like  to  achieve  the  required  accuracy 
at  the  end  of  the  "exponential"  phase.  It  is  more  advantageous  to 
overrefine  mesh  than  to  underrefine  it.  Since  the  p-version  is  using 
hierarchical  elements  [2],  as  long  as  a  mesh  is  given  the  Increase  of  N 
is  not  too  expensive,  while  for  unsmooth  solutions  a  strongly  refined  mesh 
is  very  important  get  the  desired  accuracy. 
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